The Non-Perturbative Analytical Equation of State 
for SU(3) Gluon Plasma 

V. Gogokhia and M. Vasiith 

O I HAS, CRIP, RMKI, Depart. Theor. Phys., Budapest 114, P.O.B. 49, H-1525, 

■ Hungary 

, ^ ' E-mail: gogohiaOrmki .kfki.hu, vasuthOrmki .kfki .hu 



^ 



Oh: 



o 
o 



Abstract. 

The effective potential approach for composite operators has been generahzed to 



Oj' non-zero temperatures in order to derive the analytical equation of state for pure SU{3) 



Yang-Mills fields from first principles. In the absence of external sources this is nothing 
(D ' but the vacuum energy density. The key element of this derivation is the introduction 

r^ . of a temperature dependence into the expression for the bag constant evaluated in the 

previous publications. The non-perturbative part of the analytical equation of state 

does not depend on the coupling constant, but instead introduces a dependence on the 
J^ , mass gap. This is responsible for the large-scale structure of the QCD ground state. 

r — ' Its perturbative part does analytically depend on the fine-structure constant of strong 

l/^ ■ interactions as well. As it follows from our equation of state, the two massive gluonic 

excitations with the effective masses m'^fj — 1.17 GeV and fheff = 0.585 GeV, as 
t^^ ' well as the different types of massless gluonic excitations, are present in the 5*^7(3) 



gluon plasma. Important thermodynamic quantities such as the pressure, energy and 
entropy densities, etc., have been calculated. We show explicitly that the pressure may 
continuously change around T^ = 266.5 MeV in order to achieve its Stefan-Boltzmann 
limit at high temperatures. All other thermodynamic quantities change drastically at 
/\ ' this point. The entropy and energy densities have jump discontinuities at Tc. This is 

^ . a firm evidence of the first-order phase transition in 5C/(3) pure gluon plasma. Our 

value for the latent heat is tLH — 1-54 (in dimensionless units). The heat capacity 
has a i5-type singularity (an essential discontinuity) at T^ so that the speed of light 
squared becomes zero at this point. The proposed NP analytical approach makes it 
possible to control for the first time the thermodynamics of the gluon plasma at low 
temperatures, below Tc- We have also calculated the gluon condensate, and hence the 
trace anomaly relation, as a function of temperature. Properly scaled they decrease 
not as l/T^ but as 1/T^ at high temperatures due to the explicit presence of the mass 
gap in the equation of state. All our numerical results are in very good agreement with 
corresponding lattice data at T > 2Tc. 



PACS numbers: ll.lO.Wx, 12.38.Mh, 12.38.Lg, 12.38.Aw 
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1. Introduction 

Up to now, lattice QCD remained the only practical method to investigate QCD at 
finite temperature and baryon density [H |2l [3l H]. Recently it underwent a rapid 
progress ( [H [5], El El [H |9] and references therein). However, lattice QCD, being a very 
specific regularization scheme, is primary aimed at obtaining well-defined corresponding 
expressions in order to get realistic numbers for physical quantities. One may therefore 
get numbers and curves without understanding what the physics is behind them. Such 
an understanding can only come from the dynamical theory, which is continuous QCD. 
For example, any description of the quark-gluon plasma (QGP) has to be formulated 
within the framework of a dynamical theory. The need for an analytical equation of 
state (EoS) remains, but, of course it should be essentially non-perturbative (NP), 
approaching the so-called Stefan-Boltzmann (SB) limit at very hight temperatures. 
Thus the approaches of analytic NP QCD and lattice QCD to finite-temperature QCD 
do not exclude each other; on the contrary, they should be complementary. This is 
especially true at low temperatures where the thermal QCD lattice calculations suffer 
from big uncertainties [HI HI |5l El El El [9] . On the other hand, any analytic NP approach 
has to correctly reproduce thermal QCD lattice results at high temperatures (see papers 
cited above). There already exist interesting phenomenological models based on the 
quasi-particle picture [ini[IIl[l2l[l3l[Tl[l5l[l6l[I71[T8lll9j (and references therein) to 
analyze results of S'[/(3) lattice QCD calculations for the QGP EoS. 

The main purpose of this paper is to complete the derivation of the NP analytical 
EoS for the gluon/glue plasma (GP), i.e., a system consisting purely of SU(?>) Yang- 
Mills (YM) fields without quark degrees of freedom. We have started it in our previous 
publication [20j (part I). This gluon medium was called as gluon matter (CM) in part 
I. However, here we prefer to call it GP, since we are going to include the free gluon 
contribution as well as the SB term in the present investigation, where the further 
analytical and numerical evaluation of the perturbative (PT) part of the gluon pressure 
will be also performed. The general formalism we use to generalize it to non-zero 
temperatures is the effective potential approach for composite operators [21]. In the 
absence of external sources it is nothing but the vacuum energy density (VED). This 
approach is NP from the very beginning, since it deals with the expansion of the 
corresponding skeleton vacuum loop diagrams, and thus allows one to calculate the VED 
from first principles [Sni EH [22]. The key element in this programme was the extension 
of our paper [22] to non-zero temperatures. This makes it possible to introduce the 
temperature-dependent bag constant (pressure) as a function of the mass gap. It is 
this which is responsible for the large-scale structure of the QCD ground state. The 
confining dynamics in the GP will therefore be nontrivially taken into account directly 
through the mass gap and via the temperature-dependent bag constant itself, but other 
NP effects will also be present. 

The paper is organized as follows. In sections 2 and 3 we summarize our results 
obtained in [20] for the gluon pressure at zero (i.e., the above-mentioned VED) and 
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non-zero temperatures, respectively. In section 3 the PT part of the gluon pressure is 
analytically derived as well. In section 4 we formulate a procedure how to include the 
free gluon contribution in the self-consistent way into the full GP pressure. A method of 
the simulating functions is proposed and analytical formulae for numerical simulations 
are introduced. This makes it possible to perform analytical and numerical simulations 
in order to determine the full GP pressure. In section 5 we display our numerical results 
for all the thermodynamic quantities, such as the entropy and energy densities, the heat 
capacity, etc., calculated with the help of the obtained GP pressure. We present also a 
brief discussion. In section 6 a description of the dynamical structure of the GP is given. 
In section 7 we summarize our conclusions. In appendix A the general expressions are 
given for the main thermodynamic quantities as functions of the pressure. In appendix 
B some analytical formulae for the GP thermodynamic quantities are derived, while in 
appendix C the latent heat is analytically and numerically evaluated. In appendix D 
the corresponding /3-function for the confining effective charge is explicitly derived and 
briefly discussed. 

2. The gluon pressure at zero temperature 



In order to derive the gluon pressure at zero temperature Pg in the first part of our 
investigation, we have used the effective potential approach for composite operators [21] 
to leading order, the so-called log-loop level. Analytically the gluon pressure looks like 

m 
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In the expression for the bag constant (2) symbolically shown q'^jr is the effective scale 
squared, separating the soft momenta from the hard ones in the integration over g^, 
that's < g2 < g2 



The intrinsically NP (INP) effective charge is 

.INP/„2\ 
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2 ' 
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where A^ = ^^ 

is responsible for the large-scale structure of the QCD vacuum, and thus for its INP 



AJ^y is the Jaffe-Witten (JW) mass gap ^23], mentioned above, which 
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dynamics. Let us note that how the regularized mass gap appears in QCD is shown 
in [21]. It is worth emphasizing that in our recent works [251 126] it has been proven 
that the INP effective charge (5) exactly reproduces the non-linear iteration solution 
for the full gluon propagator after the renormalization of the regularized mass gap is 
performed, so it is not an ansatz. The expression for the bag constant given in Eq. (2) 
is free of all types of the PT contributions ("contaminations"). Details of its derivation 
can be found in |22j. In the YM part of the NP gluon pressure Pym-, given in Eq. (3), 
the integration over variable q^ goes from zero to infinity. 

The PT part PpT of the gluon pressure shown in Eq. (4) contains the contribution 
which is mainly determined by the asymptotic freedom (AF) part of the PT effective 
charge, though the dependence on the INP effective charge is also present. The AF 
running effective charge a^^{q^) is given by the renormalization group equation solution, 
the so-called sum of the main PT logarithms [26l [271 EH] , namely 

^^'i<f) = . , ...!;.. .2 ^ . (6) 



l + a,61n(g2/A 



2 

YM) 



and thus like the confining effective charge (5) this is not an ansatz either. Here 
^YM = 0-09 GeV^ [29j is the asymptotic scale parameter for SU{3) YM fields, and 
b = (ll/47r) for these fields, while the strong fine structure constant is a^ = 0.1187 
[30] . Let us also note that in Eq. (6) g^ can not go below Ay^, which has been already 
symbolically shown in Eq. (4). In the g^ — > oo limit from Eq. (6) one recovers the well 
known AF expression a^'^(g^) = l/61n(g2/Ay^j) [271 EH]- The gluon pressure Pg in 
Eq. (1) is normalized to the free PT vacuum to be zero, i.e., when the interaction is 
switched formally off by putting a^^{q^) = af^^(g^) = 0, then P/vp = PpT = 0, so that 
Pg = as well. It comes from the initial normalization condition of the free PT vacuum 
to be zero in the effective potential approach up to leading order [20 

3. The gluon pressure at non-zero temperatures 



In the imaginary-time formalism [STl [32] , all the four-dimensional integrals can be easily 
generalized to non-zero temperatures T according to the prescription (let us remind that 
in part I [20] and in the present investigation the signature is Euclidean from the very 
beginning) 

l^.^T g, q' = ci' + ql = ci' + col = co' + colu;^ = 2nnT.i7) 

In other words, each integral over go of a loop momentum is to be replaced by the sum 
over the Matsubara frequencies labeled by ra, which obviously assumes the replacement 
go — 7- a;„ = 2miT for bosons (gluons). In frequency-momentum space the effective 
charges (5) and (6) become 

«f^^(g^) = «f ^(q^ ^l) = «f ^(c^^ u^l) = -^, (8) 
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and 

respectively. It is also convenient to introduce the following notations: 

T-' = /3, u = ^, (10) 

where, evidently, in all the expressions here and below q^ is the square of the three- 
dimensional loop momentum, in complete agreement with the relations (7). 

Introducing the temperature dependence into the both sides of the relation (1), we 
finally obtain 

Pg{T) = Pmp{T) + Ppt{T) = ByM{T) + Pym{T) + Ppt{T), (11) 

where P^piT) and Ppt(T) are the NP and PT parts of the gluon pressure Pg(T), 
respectively. In what follows we will call them either the NP and PT gluon pressures 
or, for simplicity, the NP and PT pressures. 

3.1. NP contribution 

One of the attractive features of the confining effective charge (8) is that it allows an 
exact summation over the Matsubara frequencies in the NP pressure Pnp{T). Collecting 
all our analytical results obtained in our previous work [20], we can write 

Pl,p(T) = ByM(T) + PYM(T) 



vr^ vr 



2 A2Pi(r) + ^T[P2(T) + P^{T) - P^{T)], (12) 



and 



while 



CO 7r^^2 f"^"^^ ^ ^ 



PiiT) = dc-— - = -T^ - / dc-,--, (13) 



^,jj e^^-l 6 Jo e^^-1' 



roo , 

P^{T) = / dco coHn (l - e-^ 

P,{T)= P"duu'ln{l-e'^-'] 

P^{T) = dojLo'^ In (1 - e'^^) , (14) 



/o 

where Uefj = 1 GeV and the mass gap A^ = 0.4564 GeV^ are fixed [201 122]. Then u' 
and uj are given by the relations 

uj' = ^q2 + 3A2 = ^oo^ + m'^ff, m'^^f = VSA = 1.17 GeV, (15) 

and 

a; = Wq2 + ^A2 = ^cu^ + rhljf, m,ff = ^A = 0.585 GeV, (16) 

respectively. The NP part of the gluon pressure Pnp{T) is shown in Fig. 1. 
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3.2. PT contribution 

In frequency-momentum space the PT part (4) of the gluon pressure (11) is 






In 1 



A + 3aJ^P{co^col] 



"^ AF/ 2 2\ 



.(17) 



However, it is instructive to begin the evaluation of the PT contribution to the gluon 
pressure (11) not from this expression, but from the very beginning, i.e., from Eq. (4). 
There is an interesting observation concerning the PT part (4). Let us consider the 
argument of its logarithm ln[l + x(g^)], where 



x{q 



4 + Sal^Piq^) (4g2 + 3A2) [1 + a^b ln(gVAf.^)] ^ 



with the help of the expressions (5) and (6). In the integral (4) g^ > Ay^, thus at its 
lower limit g^ = Ay^ the argument (18) numerically becomes 



a;(A' 



oCHs-i^YM 



YM) 



0.0185, 



(19) 



(4A|.^ + 3A2) 

and the numerical values of A^, a^ and Ay^ ^^^ ^^ given above. The argument of the 
logarithm is really small and it will become smaller and smaller with g^ going to infinity. 
This means that the logarithm ln[l + x(g^)] in the integral (4) is legitimated to expand 
in powers of a^^{q'^) or, equivalently, as. To the leading order one finally obtains 



PpT = a.A^g / 

JA 



d^q 



A^M (2vr)4 (g2 + |A2) [1 + a,Mn(gVA2.^)] 



+ 0{al). (20) 



To say that this PT term is " contaminated" by the NP contributions or to say that this 
NP term is "contaminated" by the PT contributions will be equally correct. However, 
in any case it is suppressed by the order of as in comparison with the NP part of the 
gluon pressure. In frequency-momentum space this term looks like 

q roo +°° 

PpTiT) = asA'— dojiu'T J2 

ZTX JAym »,— _o 



i^'+^ir' 



[l + asbHco^ + ul/A'yj^)] 



,(21) 



where u is given in the relation (16). For convenience, we omit the term 0{a1). 
Equivalently this integral can be presented as a sum of the two integrals, namely 



>(i) 



(2), 
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In the integral (23) the summation over the Matsubara frequencies can be performed 
analytically (i.e., exactly) with the help of formula from part I of our investigation [20] . 
Omitting all the derivation and dropping the /3-independent terms |3T], one obtains 

4')(T) = „.A^^£d..^i^, (25) 

where u is given in the relation (17). 

Thus the so-called PT contribution (22) to the gluon pressure is 

PMT) = «.A^i£ d. ." i-^ + 0(a% (26) 

and it is numerically small in the whole temperature range, indeed (see Fig. 1). 

3.3. The gluon pressure Pg{T) 

Summing up the expressions (12) and (26) up to leading order, for the gluon pressure 
Pg{T) = PNp{T) + PpT{T), (27) 

defined in Eq. (11), one finally obtains 



PAT) = ^A^ 



PiiT) + \asP[{T) 



-ryp^T) + P;{T) - P,{T)l (28) 



7r2 



where the terms of the a1 order are not shown, the integral P[{T) is 

/•oo 1 1 

PiiT) = I du uj' --^--, (29) 



U 6' 



/3d) 



and all other integrals are given in Eqs. (13) and (14). 

Let us recall that the gluon pressure (28) is normalized to be zero when the 
interaction is switched formally off, i.e., putting ctg = A^ = 0. Hence each of its 
terms in Eq. (27) becomes also zero in this case. This is due to the initial normalization 
of the free PT vacuum to zero. The confining dynamics in Eq. (28) is implemented 
via the bag pressure Bym{T) and the mass gap itself. However, other NP effects are 
also present via PyMiT) and Ppt{T), though in the PT part they are suppressed in 
the powers of as in Eq. (26), and, consequently, in Eq. (28). The gluon pressure (28) is 
shown in Fig. 1. 

4. The full GP EoS 

From our numerical results it follows that the properly scaled (i.e., divided by T^/3) 
gluon pressure (28) as a function of T has a maximum at T = 266.7 MeV, while the 
PT pressure (26) has a maximum at T = 269.4 MeV. The NP pressure (12) has 
a maximum at some finite temperature, T = Tc = 266.5 MeV. In what follows we 
will call this value " characteristic" , since the NP pressure has been calculated exactly, 
whereas the PT pressure, and hence the gluon pressure, has been calculated up to the 
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Figure 1. The gluon pressure (28), the NP pressure (12) and the PT pressure (26) 
ah properly scaled as functions of T /Tc- Effectively, all curves have a maximum at 
Tc = 266.5 MeV (vertical solid line). As functions of T they have slightly different 
maxima (see text below). 



leading order in powers of a^. Moreover, all the pressures properly scaled as functions 
of T /Tc effectively have maxima just at Tc-, see Fig. 1. From this figure it clearly follows 
that the gluon pressure (28) will never achieve the thermodynamic SB limit (A. 7) at 
high temperatures. This is not surprising, since it has been normalized to the free 
PT pressure to be zero from the very beginning [20l [211 122]. First of all, this means 
that the SB term is non-explicitly present in the gluon pressure Pg{T) (27) through the 
NP pressure P]sfp{T) (12) in the whole temperature range. From the integral P2{T) in 
Eq. (14) it can be exphcitly extracted as follows: P2{T) = J^ du cj^ln f 1 — e'^'^j = 

J^ dco u^ In (l - e-P^) - Jq^"' dco co^ In (l - e"^'^) = -(7rV8r)P5B(T) - P^(r), where 

the integral P^{T) = Jq"'' du uj'^\n(l - e~^'^). Muhiplied further by {16/n^)T, one 
finally obtains the contribution —2Psb{T) (along with terms describing another types 
of free gluons) to the NP pressure and hence to the gluon pressure. From the integrals 
P3(T) and Pi{T) in Eq. (14) it cannot be extracted explicitly in a such simple form. 
However, in the T — )■ oo limit we can neglect the difference between massless u and 
massive cu', u gluon excitations, and all the composition [P2{T) + P^iT) — P^iT)] in 
Eq. (12) goes to zero in this limit. This means that the SB pressure (A. 7) has been 
already subtracted from it, but in a very specific way. It is not surprising that the NP 
effects are important below Tc. From Fig. 1, however, it follows that they are still 
important above T,,, at least, up to 3.75Tc = 1 GeV = cUe//- The two other remarkable 
features of the gluon pressure (28) are to be also underlined. Firstly, below T^. it goes 
exponentially down and up with T going down and up, respectively. Secondly, above 
Tf. with increasing T it goes down as 1/T^ and not as 1/T^ due to the explicit presence 
of the mass gap term A^T^ (see Eq. (28) and the integral (13)), while the composition 
[P2(T) + P^iT) — P^iT)] goes to zero at high temperatures as pointed out above. In 
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general, this means that the gluon pressure (28), and hence the NP pressure (12), is 
at high temperatures a polynomial in integer powers of T up to T^. All these features 
have been discussed in more detail in |20]. The term ~ T^ in the phenomenological 
EoS, which leads to the above-mentioned fall off at high temperatures, has been first 
discussed in [33] (see also [Ml l35l 136] ). 

The gluon pressure (28) may change continuously its regime only in the near 
neighborhood of T^ in order for its full counterpart to reach the corresponding SB 
limit at high temperatures. However, the SB term cannot be simply added to Eq. (28) 
above T^, multiplied by the corresponding G-function. The problem is that the different 
types of free gluons described by the NP pressure (12) are present in the GP in 
the whole temperature range, as it was discussed above. Thus in the right-hand- 
side of Eq. (28) the free gluon contribution has to appear in a more sophisticated 
way. The most general form how this can be done is to add to Eq. (28) the term 
[e((T/T,) - l)Y{T) + e{{TjT - 1))X(T)], vahd in the whole temperature range, and 
the unknown yet functions Y{T) and X{T) are to be expressed in terms of Psb{T) and 
Pnp(T) (see subsection below). 

The previous Eq. (28) then becomes 

'^-i)x(T) + e(| 



Pgp{t) = p,(r) + e I ^ - 1 ) x(T) + e ( - - 1 ) y (T), (so) 



and its left- hand- side here and below is denoted as Pgp{T) (the above-mentioned full 
counterpart). The GP pressure (30) is continuous at Tc if and only if 

x(r,) = r(T,), (31) 

which can be easily checked. Due to the condition (31), the dependence on the 
corresponding 9-functions disappears at Tc, and the GP pressure (30) remains 
continuous at any point of its domain. The role of the function X{T) is to change 
the behavior of Pgp{T) from Pg{T) below T^ especially in its near neighborhood. The 
function Y{T) is aimed to change the behavior of Pgp{T) from Pg{T), as well as 
to introduce the SB term itself and its modification due to AF above and near Tc, 
respectively. These changes are inevitable indeed, since Pgp{T) is different from Pg{T) 
in the whole temperature range. The reason is that in the gluon pressure Pg{T) the 
SB term is missing. However, as we saw in the GP pressure it cannot be restored in 
a trivial way. That is why the appearance of the corresponding 0-functions in the GP 
pressure (30) is also inevitable, playing, nevertheless, only an auxiliary role. 

4-1. Analytical simulations 

Actual analytical and numerical simulations are the main subject of this paper in 
order to derive the NP analytical EoS for the GP valid in the whole temperature 
range. Completing this program we will be able to compare our results and achieve 
the agreement with those of the thermal QCD lattice calculations above Tc pHl [37] . 
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The space of test functions, in terms of which the functions X{T) and Y{T) should be 
found, has aheady been estabhshed. On a general ground therefore we can put 



X{T) = f,{T)PsB{T) + 0i(T)Pjvp(T), 

Y{T) = f2{T)PsB{T) + UT)P^p{T), (32) 

where all the dimensionless functions fn{T), (f)n{T), n = 1,2 will be called the simulating 
functions. The test functions Pnp(T) and Psb{T) are independent from each other, and 
they are exactly known. This also makes it possible to use in what follows the exact 
relations, derived in our previous work [20], namely 



Psb{T)-2Pnp{T)\^^^^ = Q, {A[p^^(t)-2P^p(T)]}^^^^ = 0. (33) 

In order to suppress the deep penetration of the SB term PssiT) and the additional 
NP pressure Pnp{T) below T^, the simulating functions /i(T) and 0i(T) are then 
convenient to choose as follows: 

/i(T) = Ae-"(^^/^\ a > 0, 0i(T) = Aie-"^^^^/^), ai > 0, (34) 

and A, Ai are the arbitrary constants at this stage. Evidently, these simulating functions 
mimic the asymptotic of the so-called gluon mean number Ng = [e^'^ — 1)~^ [31] in the 
/3 — !■ oo (T — )■ 0) hmit. Replacing uj by aT^ and aiTc, the different numerical factors 
A, Ai should appear in the corresponding simulating functions shown in the relation 
(34). From the GP pressure (30) below Tc then it follows 



Pgp{T) = Pg{T) + f,{T)PsB{T) + 0i(T)P^p(T) 

= PpTiT) + [1 + Aie-"i(^^/^)]P;vp(T) + Ae-"(^^/^)PsB(r), (35) 

where the relation (27) has already been used. In the T — )■ limit the additional 
contributions are exponentially suppressed, and the asymptotic of the GP pressure 
Pgp(T) is determined by the gluon pressure PgiT), as it should be. At the same time, 
these terms allow to change the value of the GP pressure from the gluon pressure near 
Tc, as it is expected. 

For the simulating function f2{T) our choice is the standard, empirical one [T0|l38]. 
recalculated only at T^ = 266.5 MeV, so it is 

fJT) = 1 - aAT) = 1 , , , , T > T,. (36) 

^^^ ' "^ ' l + 0.551n(T/Tc)' " ^ ' 

Evidently, the second term in this expression mimics the AF formula for the effective 
charge (6) as a function of temperature, while the first one leads to the correct SB limit 
for the GP pressure (30). The simulating function 02 (T) has to be a regular function of 
T as it goes to infinity; otherwise it remains arbitrary at this stage. The asymptotic of 
the GP pressure (30) at high temperatures thus becomes 
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Pgp{T) = Pg{T) + f2{T)PsB{T) + MT)Pnp{T) 



+ [1 + MT)]Pnp{T) + Ppt{T) ^ Psb{T), T^oo, (37) 

and we again use here the relation (27). Since the Ppt{T) is very small, [1 + 02(^)] has 
to be negative above T,., so the GP pressure approaches the SB limit from below, as it 
should be. The growth of the NP pressure Pnp(T) at high temperatures is suppressed 
in comparison with the growth of the SB term Psb{T) itself. The first one increases 
as A^T^ because of the explicit presence of the mass gap in the EoS (see Eq. (12) and 
integral (13)). The SB pressure increases as T^ in this limit. This means that the sum 
[1 + (l)2(T)]P]\fp(T) describes the leading and next-to-leading order term as T — )■ oo. 
We can approximate the function 4>2{T) by the corresponding expression to leading and 
next-to-leading order terms in powers of (Tc/T). 

The explicit expressions for the X{T) and Y{T) functions (32), via the chosen 
simulating functions (34) and (36), are 

X(T) = Ae-"(^=/^)P5i?(T) + Aie-"^(^^/^)Pwp(T), (38) 

and 

^'<^' - (' - l + oitV/T,) ) ^-'^) + MT)PMn (39) 

From the relations (31) and (38)-(39), and using the relations (33), it follows that 

2Ae-" + Aie-'"' = 2(1 - 0.36) + 02(Tc) = 1.28 + MTc)- (40) 

The GP pressure (30), on account of the relations (38) and (39), then looks like 

Pgp{T) = Pg{T) + f ^ - l) Ue--(^^/^)P5B(T) + Aie-"i(^^/^)p^p(T) 



T 






/-TTolf(T7^)^»(^' + ^^<^'^-<^' 



(41) 
Due the relation (40) from the previous expression at T = T^, one obtains 

Pgp{T,) = Pg(T,) + [1.28 + <P2{T,)]Pnp{T,) 

= Pg{T,) + [2Ae-" + Aie-"^]P^p(Te), (42) 

so that it depends on the number 02 (^c) only, since the values of Pnp{Tc) and Pg{Tc) 
are known from our calculations (see, for example Fig. 1). 
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4.2. Analytical and numerical simulation of (f)2{T) 

Our aim here is to find the simulating function 02 (r) by fitting lattice data at high 
temperatures above Tc f36j. For this, let us derive from the GP EoS (41) its values at 
a = T/Tc, a > 1 as follows: 

Pcpia) = BaPsBia) + Pp^a) + [1 + 02(a)]P^p(a), (43) 

where the relation (27) has been already used, and 

Ba = l- as(a) = (l ^^^^ — ) . (44) 

'^ ^ V l + 0.551nay ^ ^ 

Adjusting our parametrization of the GP pressure (43) to that used in recent lattice 
simulations for the YM SU{3) case at T = aTc in [36], one obtains 

5^ ^ ^ X 5., (45) 

where 5.26 is the general SB number shown in Eq. (A. 7) up to two digits only after the 
point. Let us note that, for simplicity, we will keep this for all our numbers through out 
this paper. So that if the calculated number, for example is 1.177..., we will approximate 
it by 1.18, if it is 1.171..., we will approximate it by 1.17, and finally if it is 1.175..., we 
will keep this value. At the same time, let us stress that we have carried out all our 
calculations keeping six digits after the point (thus we can calculate all numbers to any 
requested accuracy). Pp{T)/T'^ describes lattice data taken from the above-mentioned 
Panero's paper [SB]. They are normalized to the SB limit for the pressure (see Eq. (A. 7)), 
for example Wgp{2T^) / {2T^Y = (Pp(2Tj/(2Tj^) x 5.26 = 0.61 x 5.26 and so on. Then 
from the relation (45) and Eq. (43), one gets an equation for determining the numbers 
02(a), namely 

{p}a = 5.265, + [pt]a + {np)a[l + 02(a)], (46) 

where the numbers 

{p\a = , L, / X 5.26, {np)a = — , ',, , hot L = — , ',, (47) 

are exactly known from our calculations (see, for example Fig. 1 for the approximate 
values of the numbers {np)a and [pt\a) and 3Psb{T)/T^ = 5.26 at any high temperature 
(see Eqs. (A. 7)). The numbers {p}a are to be taken as described above. From Eq. (46) 
then it follows that 

Ma) = -1 + ^^ . (48) 

As we already know, the best way to choose the appropriate expression for the 
simulating function 02(a) is to mimic again the asymptotic of the gluon mean number 
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Ng = {e^^ — 1)^^ [31J but in the /3 — > (T — > oo) limit, which is equivalent to the 
a — 7- oo limit. Indeed, replacing u by /xTc, one can write 

A I \ -S20" B2a^ ^ _ 

so that this simulating function above T^ becomes the series in inverse powers of 
a = T/T,. 

First of all, we are interested in 

00 

02(1) = E Cn, (50) 

n=0 

as it follows from the previous Eq. (49), since 02(1) = 02 (^c)? which appears in Eq. (42). 
At values a ~ 1, we can put 

[np)a 

in Eq. (48), and the values for a in this equation are to be chosen from the low 

temperatures region, a = 1.010530, ..., 1.244705, [36]. The mean value of the numbers 

which appear in the right-hand-side of Eq. (51) is —2.12, and thus 

00 
02(1) = Ecn = -2.12. (52) 

n=0 

On the other hand, we can write 

00 m 00 m 

02(a) = E (^na-"" = E Cna-''+ E Cna-" = E Cna-^+0{a-"'^'), (53) 

n=0 n=0 n=m+l n=0 

SO that when m and a are big enough, we can approximate Eq. (53) up to leading order 

m 

02(a) = E ^na-"". (54) 

n=0 

and the values for a in this equation are to be chosen from the high temperature region, 
a = 3.006657, ..., 3.436657, [SB]. The mean value of the numbers which appear in the 
right-hand-side of Eq. (48) is —1.53, and thus 

m 

-1.53 = Ecna"". (55) 

n=0 



Then Eq. (53) can be re-written as follows: 



00 m 00 



^2(a) = Y. Cna "" = Y. (^na "" + J2 ^na "" 

n=0 n=0 n=m+l 

00 

= - 1.53 + c^+ia-"-i + E ^na-"", (56) 

n=m+2 
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which makes it possible to estimate the contribution of the remaining terms at a ~ 1. 
Indeed, in this case Eq. (56), due to the value (52), yields 



oo 



c^+i+ Y. c„ = -0.59. (57) 

n=m+2 

Substituting it back to the previous equation, one obtains 

oo oo 

02(a) = -1.53- 0.59a-'"-i-a-'"-i Y. Cn+ Y. ^nfl"''. (58) 

n=m+2 n=m+2 

This means that the whole expansion (49) or, equivalently, (58) can be effectively 
replaced by the polynomial as follows: 

02(a) = -1.53 - 0.59a-™-\ (59) 

which obviously has correct values at a = 1 and when a goes to infinity (T — )■ oo). 

The best fit to lattice data [36] above T^ is achieved at ?7i = 3 in the relation (59). 
When m < 2 the GM pressure curve lies slightly below the lattice curve in the high 
temperature limit. At m > 4 the GM pressure curve lies substantially above the lattice 
curve in the moderate temperatures region Tc < T < 2Tc, see Fig. 2. 

The final expression for the simulating function 4>2{T) is the two-term polynomial 

02(T) = -1.53- 0.59 (^^) , MTc) = -2.12, (60) 

determining this function up to leading and next-to-leading orders in the T -^ oo limit. 
Due to the value (60), from Eq. (42) one obtains 

PapiT,) = PgiT,) - 0.84P^p(T,) = 0.656 (61) 

in dimensionless units (45) and (47). At the same time, the relation (40) becomes 

v4 = -[0.42 + 0.5Aie-°^]e". (62) 

In the connection with the choice of the simulating functions a few general things 
should be made perfectly clear. In principle, all the simulating functions can be chosen 
in many different ways. However, the request that the pressure is a growing continuous 
(class C°) function of temperature in the whole range is a rather strong restriction 
(this means that it is differentiable at any point of its domain, while at T^ derivative 
itself may not be continuous, i.e., it may have a discontinuity at this point). The 
other strong condition which should be imposed on the simulating functions is that 
they should not undermine analytical and asymptotical properties of the corresponding 
test functions. The simulating functions may produce only small numerical corrections 
to the values determined by the corresponding test functions. The justified choice of 
the simulating functions can only simplify the actual numerical calculations, i.e., they 
play only an auxiliary role. Let us remind that our choice of the simulating functions 
/i(T) and 4>i{T) has been made in accordance with the asymptotic of the gluon mean 
number A''^ = (e'^'^ — 1)~^ in the /3 — )■ oo (T — )■ 0) limit. The simulating function 4'2{T) 
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Figure 2. The GP pressure (65) as a function of T/T^ is shown for different values of 
TO = 2,3,4 in the relation (59). The lattice curve [36^ for SU{2>) GP pressure is also 
shown (lower thin line). 



has been chosen in accordance with the asymptotic of the gluon mean number in the 
/3 — !• (T — )■ oo) hmit up to leading and next-to-leading orders. The gluon mean number 
describes the distribution of gluons in the GP and is closely related to its pressure and 
thus to all other thermodynamical variables. The simulating function f2{T) has been 
chosen to introduce the SB term itself and to reproduce AF corrections to it above Tc as 
well. The final result (the curve for the GP pressure (41)) will depend rather slightly on 
the numerical values of the parameters of such chosen simulating functions and details 
of their calculations. Thus this curve, and hence all the numbers behind it, will be 
uniquely determined up to the leading order (up to very small numerical uncertainties 
of the next-to- leading order). This is a general feature of the numerical simulations 
based on some analytical expressions. 

4-3. Numerical simulation of the GP pressure 

After the substitution of the relation (62), the GP pressure (41) becomes 



Pgp{T) = P,{T) 



- 
+ 6 



T 
T 



[0.42 + 0.5Aie-°^]e-"«'^=/'^)-^)P5B(T) - Aie-"^('^=/^)P;vp(T)' 
0.36 



1 



l + 0.551n(T/r, 



Psb{T) + MT)Pnp{T) 



where (t>2{T) is completely determined and given in Eq. (60). This expression is 
convenient to simulate in order to continuously (and as much smoothly as possible) 
adjust both terms at T^ which are associated with the corresponding 9-functions. 



(63) 
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It depends on the three independent (free) parameters, Ai, ai and a. The fourth 
parameter A is to be calculated through the relation (62), and thus it plays no 
independent role. The simulation procedure consists of the numerical variation of these 
parameters in order to achieve the adjustment in the requested way. In principle, there 
are no other constraints on these constants apart from that all the derivatives of the GP 
pressure (63) should not gain negative values at low temperatures below T^, i.e., they 
should exponentially approach zero from above. 
Our best values for them are: 

Ai = -1.1, ai = 0.006, a = 4. (64) 

The value a = 4 means that the deep penetration of the SB term below T^ is indeed 
suppressed, as expected. The value ai = 0.006 means that the NP pressure is not 
suppressed, but only the slope of its curve is decreasing below T^, which is also expected. 
The number Ai is, in principle, a dimensionless combination of all our dimensional 
parameters: A^, oOeff, T^, AyM- Let us remind that only A^ or, equivalently, We// and 
Ay^j are independent, while the latter one is suppressed by a^, see Eq. (26). However, 
how to assign to it a physical meaning (if any) is not clear; anyway this is not important. 
Thus the GP pressure (63) finally becomes 



0.146e-^«^=/^)-i)P5B(T) - l.le-°-°°6(^^/^)P;vp(T) 



Pgp{T) = P,{T) 




-(^ 


-1 


-(I- 


-1 



/- i + oitV/rJ ^--'^'-^^-'^'^--'^' 



(65) 



and it is completely determined now, since 4'2{T) is also determined and given by 
the relation (60). All other thermodynamic quantities will be calculated just via this 
expression. 

5. Numerical results and discussion 

In Fig. 3 we present the GP pressure (65) as well as the lattice pressure [36]. The 
main difference between them lies in the region of low temperatures below Tc. Within 
the NP analytic approach the confining effective charge (8) is mainly responsible for the 
structure of the GP in this region. It gives rise to the gluonic excitations of the dynamical 
origin which cannot be taken into account by the lattice simulations. The corresponding 
/3-function (see appendix D) is always negative, as it is required by confinement [26|l27]. 
At the same time, the lattice /3- function is always positive and almost constant [361 EZ] • 
The discrepancy between NP analytic GP and lattice pressures seen in Fig. 3 just above 
Tc up to 2Tc is also due to other NP effects which are still important in this region as 
pointed out above (see discussion in connection with Fig. 1). 
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Figure 3. The GP pressure (65) as a function of T/Tc (upper tliick line). The lattice 
data [36] for SU{'ii) GP pressure is also shown (lower thin line). Both pressures are 
scaled in the same way in accordance with Eq. (45). 



The GP entropy and energy densities are shown in Fig. 4. 
discontinuity in the energy density, the so-called latent heat (LH) is 

eLH = 1.54 



The size of the 



(66) 



in dimensionless units (see appendices B and C for its definition and analyti- 
cal/numerical derivation, respectively). This means that the first-order phase transition 
in the GP is analytically confirmed for the first time, in complete agreement with ther- 
mal 5*^7(3) YM lattice simulations. The physical reason of such sharp changes at T^ in 
the derivatives of the GP pressure is that its exponential rise below Tc is changing to 
the polynomial rise above Tc in order to reach finally the thermodynamic SB limit. Our 
value (66) for the LH is consistent with various thermal SU{3) YM lattice calculations: 
in [37] eLH ~ 2, in [36] eiH = 1-4, in [39] eLH = 1-413 and in [ID] eLH = 1-54. Note, 
that our value is slightly bigger than the recent lattice values in [361 ES] ) and coincides 
with the value in [lO] (for simplicity, we do not specify the lattice parameters such as 
Nt, Ng, etc.). This agreement is rather remarkable, taking into account the completely 
different natures of the NP analytical and lattice approaches to thermal QCD. 

The trace anomaly relation, which is especially sensitive to the NP effects including 
the conformal symmetry breaking, is shown in Fig. 5. Its two most interesting features 
are: the rapid rise of the peak due to the LH in the energy density (see Fig. 4) is exactly 
placed at T^, while in lattice calculations it peaks at about l.lTc [351 EHl ISH HO]. We 
interpret this lattice artefact as a clear evidence that by lattice simulations one cannot 
get correct numbers for all the thermodynamic quantities near Tc and at T^ itself. As 
it follows from Fig. 3, the lattice pressure below T^ is almost zero and at T^ it is about 
0.1, while our pressure is 0.656 (in dimensionless units). That was the reason why we 
simulate the function 4>2{T) in subsection (4.2) not using the lattice value at Tc. It is 
worth underlining once more that lattice simulations cannot take into account the NP 
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Figure 4. The GP pressure, the entropy and energy densities all properly scaled are 
shown as functions of T/Tc- The finite jumps in densities are clearly seen, and the LH 
is ^LH ~ 1.54. Their common SB limit (A. 7) is approaching rather slowly. 




Figure 5. The trace anomaly properly scaled is shown as a function of T/Tc (solid 
line) . The gluon condensate at zero temperature scaled in the same way is also shown 
(dot-dashed line). 



dynamical effects, especially confinenient, neither at zero nor at non-zero temperatures. 
They dominate the structure of the QCD vacuum at zero temperature and will remain 
important at non-zero temperatures as well (up to 3.75Tc within our approach [20]). In 
other words, they cannot be so negligible small below Tc as it is shown by all lattice 
results for all the thermodynamic quantities. However, above T^. at sufficiently high 
temperatures T > 2Tc our results coincide because of the adjusted procedure necessary 
performed in subsections 4.2 and 4.3. The second interesting feature is that like in all 
lattice simulations (see for example [351 ESI |371 HO]) the trace anomaly falls off as 1/T^ 
and not as 1/T^ at high temperatures. As emphasized above in section 4 of the present 
investigation and in [20] , this is due to the explicit presence of the mass gap term A^T^ 
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Figure 6. The properly scaled gluon condensate at non-zero temperatures is shown 
as a function of T/Tc. It approaches zero from below at high temperatures. 




Figure 7. The NP gluon condensate scaled by T^ as a function of T jT^. The zero 
temperature gluon condensate scaled by T** is also shown. At high temperatures the 
NP gluon condensate scales as 1/T^. 



in the EoS (65). The terms proportional to T^ in the trace anomaly cancel each other 
at high temperatures due to Eq. (A.9), and hence the mass gap term comes out into 
play in this limit. It is the mass gap which violates the conformal symmetry from the 
very beginning within our approach. 

In close connection with the trace anomaly is the gluon condensate defined by 
Eq. (A. 6) and shown in Fig. 6. It approaches zero from below, so it gains negative values 
at high temperatures (fixed also by lattice simulations in [37]) due to the behavior of the 
trace anomaly in this region. In this connection, a few remarks are in order. The gluon 
condensate at zero temperature has been defined and calculated by the subtraction of 
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Figure 8. P(e) EoS (solid line) and approach to con formality = (1/3) (thin line). 




Figure 9. Conformity C{T) = P{T)/e{T) as a function of T/Tc. It is zero at T = 0, 
and has a jump at Tc due to a jump in the energy density at this point. It shows 
the non-trivial dependence on T below Tc, and approaches the SB limit 1/3 at high 
temperatures rather quickly. 



all types of the PT contributions ("contaminations") in ^22j. Here the temperature- 
dependent gluon condensate has been calculated with the help of the GP pressure (65), 
which contains explicitly the PT SB terms in different combinations. If the temperature- 
dependent gluon condensate is to be defined in the same way as its zero-temperature 
counterpart, then only the NP pressure (12) should be used for its calculation. This 
has been already done in [20] and shown in Fig. 7. It is always positive and at high 
temperatures it scales as 1/T^. 

The GP pressure versus the GP energy density, the so-called EoS -P(e), is present 
in Fig. 8. The size of the LH and the rather quick approach to conformality are clearly 
seen. Let us note in advance that we distinguish between the conformality here and 
conformity defined in Eq. (A. 4), though numerically in the limit of high temperatures 
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Figure 10. The heat capacity as a function oiT/Tc- It has the (5-type singularity (an 
essential discontinuity) at Tc- It approaches its SB limit (A. 7) (horizontal dashed line) 
at high temperatures very slowly. 




Figure 11. The velocity of sound squared as a function of T/Tc. It shows the non- 
trivial dependence on T below Tc, while at T = and aXT = Tc it is zero. It approaches 
the SB limit 1/3 at high temperatures rather quickly. 



they are the same. The conformity itself is shown in Fig. 9. It has a finite jump at 
Tc because of a jump in the energy density at this point, and approaches its SB limit 
(A. 8) at high temperatures rather quickly. However, its most interesting feature is its 
non-trivial dependence on T below T^, which has been fixed explicitly in SU{?)) GP for 
the first time. It is due to the fact that the conformity is the ratio of the independent 
thermodynamic quantities, and therefore is not obliged to be always a growing function 
of T below Tci though above Tc it is always growing. 

The last independent thermodynamic quantity, the heat capacity defined in 
Eq. (A. 2), is shown in Fig. 10. It is always a growing function of T below and above T^ 
while at Tc it has a 5-type singularity (an essential discontinuity) due to the expression 
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(B.6). It approaches its SB limit (A. 7) at high temperatures very slowly. The velocity 
of sound squared is shown in Fig. 11. Like conformity it is the ratio of the independent 
thermodynamic variables, see Eq. (A. 3). It behaves very similar to conformity (Fig. 9), 
since the latter one mimics its properties. The principal difference from conformity is 
that at Tc it is zero because of the heat capacity having the above-mentioned 5-type 
singularity at this point. It approaches its SB limit (A. 7) rather quickly. 

6. The dynamical structure of SU(?)) GP. A brief description 

A rich dynamical structure of the GP emerges in the framework of our approach. The 
NP part (12) of the GP pressure (65) tells us that in the whole temperature range 
we have two different massive gluonic excitations (or, equivalently, effective gluonic 
degrees of freedom ITOl ST]) u' and u with the effective masses Tn'^^j = 1.17 GeV and 
ffieff = 0.585 GeV, respectively. Both effective masses are due to the mass gap A^, which 
is responsible for the large-scale structure of the QCD ground state. It is dynamically 
generated by the nonlinear interaction of massless gluon modes [211 |25l |26]. The first 
massive excitations can be interpreted as glueballs, since rn'^jf is compared to the masses 
of scalar glueballs [12]. The second one frieff, might be identified with an effective gluon 
mass of about (500 — 800) MeV, which arises in different quasi-particle models (see again 
|42] and references therein). 

We also have two different massless gluonic excitations u: the NP massless gluons, 
propagating in accordance with the integral -Pi(T) in Eq. (13), and almost SB massless 
gluons, propagating in accordance with the integral P2{T) in Eqs. (14). Let us remind 
that the integral (13) should be multiplied by (6/'7r^)A^, and all other integrals (14) are 
to be multiplied by {16/tt^)T, when one speaks about different NP contributions to the 
pressure (12). However, this is not the whole story yet. The NP pressure (12) enters the 
GP pressure (65) in three different places, one through the gluon pressure (27). This 
means that all the above-mentioned massive and massless gluonic excitations propagate 
in the GP in two different ways: below and above T^ they propagate in accordance with 
the terms in Eqs. (35) and (37) with parameters and the function 4'2{T) fixed by the 
relations (64) and (60), respectively. 

The PT part (26) of the gluon pressure (27) assumes the presence of the NP 
massive gluonic excitations u in the GP in the whole temperature range. Having the 
same effective masses, nevertheless, they propagate in a completely different way, which 
is determined by the mass gap A^ and Aym (compare the integrals (26) and (14)). 
However, their propagation is suppressed by the order of a^. The terms of the order 
of al and higher have been left undetermined. In any case, numerically they are very 
small. Let us emphasize once more that all these NP massive and massless gluonic 
excitations have not been introduce by hand; contrary, they are of the dynamical origin, 
and only they (i.e., of the dynamical origin) can be described and accounted for within 
our approach to thermal QCD. 

In the GP in the whole temperature range the free gluons context within our 
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approach is to be described as follows: we definitely have the SB free massless gluons 
far away from T^, while near T^. their propagation is corrected by AF. Below T^. the deep 
penetration of the SB term itself is suppressed. So we have three different types of the 
PT free gluons, apart from those discussed above and at the beginning of section 4. 

The exponential rise of all the independent thermodynamic quantities in the 
transition region (0.6 — l)Tc [20] clearly seen in Figs. 4 and 10 indicates that near 
Tc a dramatic increase in the number of effective gluonic degrees of freedom will appear 
(for example, the glueballs will begin to dissolve rapidly). This will lead to drastic 
changes in the structure of the GP. A change in this number is enough to generate 
pressure gradients, but not enough to affect the pressure itself. It varies slowly and 
therefore remains continuous in this region. At the same time, the pressure gradients 
such as the energy and entropy densities, etc., undergo sharp changes in their behavior, 
having different types of discontinuities at Tc. Thus SU{?>) GP has a first-order phase 
transition within our approach, in complete agreement with the thermal QCD lattice 
calculations (see [361 [37] and references therein). Of course, not all the glueballs will be 
dissolved in the transition region. Some of them will remain above T^, along with other 
massive and massless gluonic excitations, forming thus a possible mixed and extended 
mixed phases around T^ from its both sides [20| 113] . One can conclude that the NP 
physics of the mixed phases (the temperature interval approximately (06 — 3.75)Tc) in 
S'f/(3) GP is well understood within our approach. At very high temperatures starting 
with 3.75T'c the NP effects become very small, and the structure of the GP will be mainly 
determined by the SB relations (A.7)-(A.9) between all the thermodynamic quantities. 

7. Conclusions 

The effective potential approach for composite operators [21] has been generalized to 
non-zero temperatures in order to derive the NP analytical EoS for pure S'[/(3) YM fields 
from first principles. In its NP part (12) there is no dependence on the coupling constant, 
only a dependence on the mass gap, which is responsible for the large-scale structure of 
the QCD ground state. A key element of this part is the generalization of the expression 
for the bag constant at zero temperature [201 122] to non-zero temperatures. Its PT part 
does analytically depend on the fine-structure constant of strong interactions as well. 
We also propose and formulate a method how to include the free gluon contribution into 
the GP pressure (65) in a self-consistent way. 

Our main quantitative and qualitative results in this investigation are: 

(i). In the GP EoS (65) the confining dynamics at non-zero temperatures (8) is 
taken into account through the T-dependent bag constant and the mass gap A^ itself. 
The GP pressure (65) compared to lattice pressure [36] is shown in Fig. 3. 

(ii). Other NP effects are also taken into account via the YM and PT parts of the 
GP pressure, though the latter one is suppressed by the order of a^. 

(ill). The mass gap A^ or, equivalently, cUe// and the asymptotic scale parameter 
for YM fields Ay^j are only two independent input scale parameters in our approach. 
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(iv). The characteristic temperature T^ = 266.5 MeV is a temperature at which 
the maximum of the NP part of the GP pressure is achieved (Fig. 1). 

(v). The presence of the mass gap A^T^ term in the GP pressure (65). 

(vi). Due to this the trace anomaly and the NP gluon condensate will go down as 
1/T^ at high temperatures. The trace anomaly has a peak just at Tc (Figs. 5,6 and 7). 

(vii). The presence of the two different types of massive gluonic excitations of 
the NP origin with an effective masses 0.585 GeV and 1.17 GeV. They propagate in 
different ways below and above T^. There is no place for some other types of massive 
gluonic excitations in our picture. 

( viii) . The presence of the different types of massless gluonic excitations including 
the free PT or, equivalently, SB gluons far away from T^. 

(ix). Below Tc all the independent thermodynamic quantities exponentially go 
down and up with T going down and up, respectively (Figs. 4 and 10). Above Tc this 
regime is changed to a polynomial growing one. 

(x). The thermodynamic quantities being the ratios of the corresponding 
independent ones such as conformity and the velocity of sound squared demonstrate the 
non-trivial similar dependence on T below Tc (Figs. 9 and 11). At high temperatures 
their approach to the SB limit (A. 8) is accelerated by a factor of 3. 

(xi). At Tc the GP pressure is continuous, while the energy and entropy densities 
have finite jumps. The size of the discontinuity in the energy density is: eiH = 1-54 (in 
dimensionless units), see Fig. 4. The heat capacity has an essential discontinuity of a 
S-tjpe singularity at Tc, see Fig. 10. 

(xii). The first-order phase transition at Tc in SU{3) GP is thus analytically proven. 

(xiii). Since the NP effects are still important up to 3.75Tc, all the independent 
thermodynamic quantities approach their perspective SB limits rather slowly. 

(xiv). Because of this the behavior of SU{3) GP in this region is substantially 
different from the behavior of a gas of free massless gluons. 

The NP analytical approach to thermal QCD provides a detailed numerical and 
dynamical description of SU{3) GP in the whole temperature range, and especially 
in the region of low temperatures below Tc. Something that's not accessible for and 
missing by lattice approach to thermal QCD. Our next and final step is to include 
the quark degrees of freedom into the formalism of the effective potential approach for 
composite operators at non-zero temperatures. This will make it possible to derive the 
NP analytical EoS for the QGP. 
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Appendix A. Main thermodynamic quantities 

Together with the pressure P{T), the main thermodynamic quantities are the entropy 
density s{T) and the energy density e(T). The general formulae which connect them 
are [31] 

e(T) = T (^^) - P{T) = Ts{T) - P{T) (A.l) 

for pure YM fields, i.e., when the chemical potential is equal to zero. Let us note that 
in quantum statistics the pressure P{T) is nothing but the thermodynamic potential 
n{T) apart from the sign, i.e., P{T) = -fi(r) > 0. 

Other thermodynamic quantities of interest are the heat capacity cv{T) and the 
velocity of sound cj.(T), which are defined as follows: 

and 

^,^T) -'^-^y (A.3) 

de{T} cv[T) 

i.e., they are defined through the second derivative of the pressure. The so-called 
conformity 

C(T) = ^ (A.4) 

mimics the behavior of the velocity of sound squared (A.3) but without involving such 
a differentiation. 

A thermodynamic quantity of special interest is the thermal expectation value of 
the trace of the energy momentum tensor. This trace anomaly relation measures the 
deviation of the difference 

e(T) - 3P(T) (A.5) 

from zero at finite temperatures, in the high temperature limit it must vanish. As a 
consequence it is very sensitive to the NP contributions to the EoS. It also assists in the 
temperature dependence of the gluon condensate [44j (see [37] as well), namely 

< G^ >T=< G^ >o -[e(T) - 3P(T)], (A.6) 

where < G^ >o=< G"^ >t=o denotes the gluon condensate at zero temperature, whose 
numerical value is discussed in appendix B of [20] . 
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The high-temperature behavior of all the thermodynamic quantities is governed by the 
SB ideal gas limit, when the matter can be described in terms of non- interacting massless 
particles (gluons). In this limit these quantities satisfy special relations [31] 

WsBiT) _ tsBiT) _ 3ssBiT) _ cvisB)iT) _ 24 _2 



ji4 



J^4 



47-3 



4T3 



45 



and 



Csb{T) = cl^SB)iT) = ^, T^oo, 



vr^ ^ 5.26, T -^oo, (A.7) 



(A.8) 



on account of the previous relations and their definitions in Eqs. (A. 3) and (A. 4). The 
trace anomaly relation (A. 5) also satisfies the SB limit, namely 



esB{T)-3PsB{T) = 0, T ^ oo. 



(A.9) 



as it comes out from the relations (A.7). In what follows its right-hand-side will be 
called the general SB number. 

Appendix B. Analytical formulae for the GP thermodynamic quantities 

It is instructive to derive analytically all the necessary formulae for the thermodynamic 
quantities using the GP pressure (30). Differentiating it in accordance with the definition 
(A.l), on account of the relation (31), one obtains 



s{T) 



dT 



+ 0(^ 



dX{T) 
dT 



+«i|-^ 



dY{T) 
dT ' 



(B.l) 



for convenience, the subscript "GP" in all the thermodynamic quantities will be omitted 
here and below. It is easy to see that the entropy density has a jump at Tc, 



As{T,) = [s{T > T,) - s(r < TJ] 



T~^Tc 



dY{T) dX{T) 



dT 



dT 



(B.2) 



Jr=r, 



where the difference in the right-hand-side of this equation has to be positive. 
In the same way for the energy density, one obtains 



dP (T) 

<T) = T— 1^ - p,{T) + e 



dT 
+ 0(^-1 



T 



T^-A-m 



r'-^-rm 



The size of the discontinuity in the energy density (latent heat (LH)) is 

'dY{T) 



eLH{T,) = Ae{T,) = [e{T > TJ - e(T < T,)] 



T^Tc 



T 



(B.3) 



dX{T) 



dT 



dT 



[B.4) 



Jt=Tc 
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and thus it is in agreement with the discontinuity in the entropy density, since from 
Eqs. (B.3) and (B.5) it follows that 



eLH{Tc) = ^e{T,) =TAs{T,). 



(B.5) 



The next independent thermodynamic quantity is the heat capacity defined in 
Eq. (A. 2). Differentiating the entropy density (B.l), one finally obtains 

The important observation is that the heat capacity has a (5-type singularity (an essential 
discontinuity) at T = T^, so that the speed of light squared (A. 3) at this point is zero 



cliZ 



cv(Tc 



0, 



(B.7) 



indeed. All other thermodynamic quantities can be analytically derived by using these 
formulae and the corresponding definitions given in the previous appendix A. 

Appendix C. Analytical and numerical evaluation of the latent heat 

From the expression (38) it follows that 



^ = «A||e-"(-=/-)p,,(T) + Ae-"(-=/-)^^ 



+ aiAi||e--(^^/^)p^p(T) + A,e^^^(T.mdP^^ 



and at T = Tf. it is 



(C.l) 






[-^)^-^^'^PssiT.) + Ae 



\ dT 



Tc 



+ ^e-P.p(T.) + A,e-f^^-^' 



T. 



dT 



(C.2) 



In the same way from the expressions (39) and (60) it follows that 



dY{T) 
dT 



0-36 0.55 .^., /,_ 

:i + 0.551n(T/Tj)2 j. sbk ) ^ y ^ 



0.36 



0.551n(T/Tc 



+ 



2.36 fT, 



dPNpiT) 



^^ PMT) + MT)- ^ 



dPsBJT) 
dT 

(C.3) 



and at T = Tc it is 
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2.36 ,^, , ^,^J dPr,p{T) \ 
+ -j^Pnp{Tc) + MTc) I — ^ 1 • (C.4) 

Taking into account these relations and the relations (33) and (40), the latent heat 
(B.4) thus becomes 

eLH{Tc) = [2.756 - 2aAe-" - a,A,e~'''] Pnp{T,), (C.5) 

and in dimensionless units eiu = {^lh(Tc)/T^) it is 



^LH 



2.756 - 2aAe-" - aiAic'"' x 0.875, (C.6) 



where the number Pi^p{Tc)/T^ has been already substituted. Using further the numbers 
(64), one finally obtains 

eLH = 1.54. (C.7) 

Appendix D. The /3-function for the confining effective charge at non-zero 
temperatures 

It is instructive to show explicitly the corresponding /3-function for the INP effective 
charge (5). From the renormalization group equation, 

r ' ^g'^^ =/3(aing^;A^)), (D.l) 

it simply follows that 

Thus, the corresponding /3-function as a function of its argument is always in the domain 
of attraction (i.e., negative). So it has no IR stable fixed point indeed as it is required 
for the confining theory [2^. In frequency-momentum space from Eqs. (8) and (D.2) 
one gets 

nINP(,2 , 2\ ^INP(,2 , 2\ ^ /'■Pl Q^ 



/3(af ^(g^; A^)) = -af ^(g^; A^) = -^. (D.2) 



00^ + UJ. 



n 



The confining effective charge with the corresponding /3-function (D.3) determines the 
structure of SU{3) GP at low temperatures as well as in the extended mixed phase up 
to 3. 75 Tc [201 SS] within our approach. 
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